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is ABSTRACT 

19 We report on the VERITAS discovery of very-high-energy (VHE) gamma- 

ray emission above 200 GeV from the high-frequency-peaked BL Lac object 
RXJ0648. 7+1516 (GB J0648+1516), associated with 1FGL J0648. 8+1516. The 
photon spectrum above 200 GeV is fit by a power law dN/dE = F (E / 'E )~ r 
with a photon index T of 4.4 ± 0.8 sta t ± 0.3 syst and a flux normalization F of 
(2.3 ±0.5 stat ± 1.2 sys ) x 10- 11 TeV^cm-V 1 with E = 300 GeV. No VHE vari- 
ability is detected during VERITAS observations of RX J0648. 7+1516 between 
2010 March 4 and April 15. Following the VHE discovery, the optical identifica- 
tion and spectroscopic redshift were obtained using the Shane 3-m Telescope at 
the Lick Observatory, showing the unidentified object to be a BL Lac type with a 
redshift of z = 0.179. Broadband multiwavelength observations contemporaneous 
with the VERITAS exposure period can be used to sub-classify the blazar as a 
high-frequency-peaked BL Lac (HBL) object, including data from the MDM ob- 
servatory, Swift-WOT and XRT, and continuous monitoring at photon energies 
above 1 GeV from the Fermi Large Area Telescope (LAT). We find that in the 
absence of undetected, high-energy rapid variability, the one-zone synchrotron 
self-Compton model (SSC) overproduces the high-energy gamma-ray emission 
measured by the Fermi-LAT over 2.3 years. The SED can be parameterized sat- 
isfactorily with an external- Compton or lepto-hadronic model, which have two 
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and six additional free parameters, respectively, compared to the one-zone SSC 
model. 

20 Subject headings: gamma rays: galaxies - - BL Lacertae objects: individual 

(RXJ0648.7+1516, 1FGL J0648. 8+1516, VER J0648+152) 



1. Introduction 



23 1FGL J0648. 8+1516 was detected b y Fermi-LAT in the first 11 months of operation at 

24 greater than 10 standard deviations, a (lAbdo et al.ll2010a[ ). This source was flagged as a 

25 very-high-energy (VHE; E>100 GeV) emitting candidate by the Fermi-LAT collaboration 

26 by searching for >30 GeV photons. This information triggered the VERITAS observations 

27 reported here. 1FGL J0648.8+151 6 is found to be assoc iated with RXJ0648. 7+1516, which 
2B was first discovered by ROSAT (IBrinkmann et al.lll997l ). A radio counterpart was identified 
29 in the NRAO Green Bank survey (IBecker et al.lll99ll). Two s ubsequent attempts to identify 



30 an optical counterpart were unsuccessful (IMotch et al.l 11998c lHaakonsen et al. 



20091 ). 



At 6° off the Galactic plane and without optical spectroscopy, the nature of this object 
remained unknown until optical spectroscopy was obtained in response to the VERITAS 
detection. These observations allow the active galactic nucleus (AGN) to be classified as a 
BL Lac, a type of AGN that has a jet co-aligned closely with the Earth's line of sight and 
displays weak emission lines. These AGN are characterized by non-thermal, double-peaked 
broadband spectral energy distributions (SED). Based on the radio an d X-ray flux, the BL 
Lac c an further be classified as a high-frequency-peaked BL Lac (HBL) (IPadovani fc Giommi 



19951 ). or if classified by the lo cation of its low-energy peak, a high-synchrotron-peaked BL 



39 Lac (HSP) flAbdo et al.ll2010bh . 



2. Observations and Analysis 



2.1. VERITAS 



VERITAS comprises four imaging atmospheri c Cherenkov telescopes and is sensitiv e 
to gamma-rays between ~100 GeV and ~30 TeV jWeekes et all I2OO2I : rHolder et alJboOfih . 
The VERITAS observations of RX J0648. 7+1516 were completed between 2010 March 4 and 
April 15 (MJD 55259-55301), resulting in 19.3 hours of quality-selected live time. These 
observations were taken at 0.5° offset in each of four directions to enable simultaneous back- 



ground estimation using the reflected-region method (jFomin et al.lll994l ). 



- 5 - 



48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 

61 
62 
63 
64 
65 

66 
67 
68 
69 
70 
71 



The VERITAS events are parameterized by the principal moments of the elliptical 
shower images, allowing cosmic-ray background rejection through a set of selection criteria 
(cuts) which have been optimized a priori on a simulated, soft-spectrum (photon index 4.0) 
source with a VHE flux 6.6% of that observed from the Crab Nebula. The cuts discard 
images with fewer than ~50 photoelectrons. Events with at least two telescope images 
remaining are then cosmic-ray discriminated based on the mean-scaled-width (MSW) and 
the mean-scaled-length (MSL) parameters. Events with MSW < 1.1, MSL < 1.4, a height of 
maximum Cherenkov emission > 8 km and an angular distance to the reconstructed source 



( Coean 


2008; 


Daniel 


2008) 



After background rejection, 2711 events remain in the source region, with 16722 events 
remaining in the background regions (larger by a factor of 6.89) . The 283 exc ess events 



result in a significance of 5.2<r, calculated using Equation 17 from lLi fc Mai (119831 ). 



A differential power law dN/dE = F (E /300 GeV)" r is fit to the VERITAS data 
from 200 to 650 GeV, shown in the top panel of Figure 1. The fit (x 2 = 0.90 with 3 
degrees of freedom (DOF), probability of 0.83) results in a flux normalization of F Q = (2.3 ± 



0.5 stai ± 1.2 syst ) x 10 



-ii 



photons cm 2 s 1 TeV 1 and an index of T = 4.4 ± 0.8 sta t ± 0.3 



syst) 



corresponding to 3.3% of the Crab Nebula flux above 200 GeV. 

The angular distribution of the excess events is consistent with a point source now 
designated VER J0648+152, located at 102.19° ± 0.11° tat RA and 15.27° ± 0.12° tat Dec 
(J2000). The systematic pointing uncertainty of VERITAS is less than 2 5" (7xl0~ 3 degrees) . 



This position is consistent with the radio position of RXJ0648. 7+1516 (IBecker et al.lll99ll ). 
A nightly-binned VHE light curve is fit with a constant and shows a \ 2 nu U hypothesis 
probability of 0.39, showing no significant variability during the observation. 



2.2. Fermi-L AT 



73 The Fermi-LAT is a pair- conversion teles c ope sensitive to p hotons between 20 MeV 

74 and several hundred GeV ( Atwood et al. 2009 ; Abdo et al. I bood l The data used in this 

75 paper encompass the time interval 2008 Aug 5 through 2010 Nov 17 (MJD 54683-55517), 

76 and were analyzed with the LAT ScienceTools software package version v9rl5p6, which 

77 is available from the Fermi Science Support Center (FSSC). Only events from the "diffuse" 

78 class with energy above 1 GeV within a 5° radius of RX J0648. 7+1516 and with a zenith 

79 angle < 105° were used. The background was parameterized with the files gll_iem_v02.fit and 
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so isotropic_iem_v02.txt 0. The normalizations of the components were allowed to vary freely 

si during the spectral point fitting, which was performed with the unbinned likelihood method 

82 and using the instrument response function P6_V3_DIFFUSE. 

83 The spectral fits using energies above 1 GeV are less sensitive to possible contamination 

84 from unaccounted (transient) neighboring sources, and hence have smaller systematic errors, 

85 at the expense of slightly reducing the number of source photons. Additionally, there is no 
se significant signal from RXJ0648. 7+1516 below 1 GeV. The analysis of 2.3 years between 
87 2008 Aug 5 and 2010 Nov 17 (MJD 54683-55517) of Fermi-LAT events with energy between 
sa 0.3-1 GeV (fixing the spectral index to 1.89) yields a test statistic (TS) of 9, corresponding 

89 to ~ 3a B In addition to the background, the emission model includes two nearby sources 

90 from the 1FGL catalog: the pulsars PSR J0659+1414 and PSR J0633+1746. The spectra 

91 from the pulsars are parameterized with power-law functions with exponential cutoffs, and 

92 the values are fixed to the values found from 18 months of data. The spectral fluxes are 

93 determined using an unbinned maximum likelihood method. The flux systematic uncertainty 

94 is estimated as 5% at 560 MeV and 20% at 10 GeV and abovefl 

95 The results from the Fermi-LAT spectral analysis are shown in the bottom panel of 

96 Figure 1. There is no variability detected in four time bins evenly spread over the 2.3 

97 years of data. The dataset corresponding in time to the VERITAS observations between 

98 between 2010 March 4 and April 15 (i.e. MJD 55259—55301) does not show any significant 

99 signal and thu s we report 2a upper limits that were computed using the Bayesian method 



100 ( iHelene Ill983[ ). where the likelihood is integrated from zero up to the flux that encompasses 

101 95% of the posterior probability. When using the data accumulated over the expanded full 

102 2.3 years of data, we find that 1FGL J0648. 8+1516 is significantly detected above 1 GeV 

103 with a TS of 307. The spectrum is fit using a single power-law function with photon flux 

104 FyiGeV — (1-8 =t 0.2 s t a t) x 10~ 9 photons cm _2 s _1 and hard differential photon spectral index 

105 r lat = 1-89 ± 0.10 stat . The analysis is also performed on five energy ranges equally spaced 
we on a log scale with the photon index fixed to 1.89 and only fitting the normalization. The 

107 source is detected significantly (TS>25) in each energy bin except for the highest energy 

108 (100-300 GeV), for which a 95% confidence level upper limit is calculated. 



1 The files are available at http : //f ermi . gsfc .nasa.gov/ssc/data/access/lat/BackgroundModels .html 

2 See lMattox et al.1 (jl996T ) for TS definition. 



See http : //f ermi . gsfc . nasa. gov/ ssc/data/analysis/LAT_caveats .html 
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2.3. Swi/i-XRT 



no The Swift-XRT (IGehrels et al.ll2004[ iBurrows et al.ll2005l ) data are analyzed with HEA- 

iii SOFT 6.9 and XSPEC version 12.6.0. Observations were taken in photon counting mode 

n2 with an average count rate of ~ 0.3 counts per second and did not suffer from pile-up. Six 

m target-of-opportunity observations summing to 10.5 ks were collected on six different days 

H4 between 2010 March 18 and April 18 (MJD 55273 and 55304), inclusive. These observations 

us were combined with a response file created from summing each observation's exposure file 

lie using ximage. The photons are grouped by energy to require a minimum of 30 counts per 

H7 bin, and fit with an absorbed power law between 0.3 and 10 keV, allowing the neutral hy- 

n8 drogen (HI) column density to vary. A HI column density of 1.94 ±0.14 x 10 21 cm ~ 2 is found, 

n9 only slightly higher than the 1.56 x 10 21 cm -2 quoted in Kalberla et al.l ( 120051 ). The com- 

120 bined X-ray energy spectrum is extracted with a fit (x 2 = 114 for 88 DOF, null hypothesis 

121 probability of 3.2xl0~ 2 ) with a photon index of 2.51 ±0.06 and an integral flux between 0.3 

122 and 10 keV of (1.24 ± 0.03 stat ) x 10~ n ergs cm" 2 s -1 . This corresponds to a 0.3 to 10 keV 

123 rest frame luminosity of 1.1 x 10 45 ergs s" 1 . The deabsorbed spectrum is used to constrain 

124 modeling. 



2.4. Swift-UVOT 



126 The Swift-XRT observations were supplemented with UVOT exposures taken in the U, 

127 UVM2, and UVW2 bands (cen tered at 8.56 x 10 14 Hz, 1.34 x 10 15 Hz, and 1.48 x 10 15 Hz, 
i2B respectively; iPoole et al.l ( 120081 )). The UVOT photometry is performed using the HEASOFT 

129 program uvotsource. The circular source region has a 5" radius and the background regions 

130 consist of several circles with radii between 10 — 15" of nearby empty sky. The results are 

131 reddening corrected using R(V)=3.32 and E(B- V)=0.14 (jSchlegel et al.lll998f) . The Galactic 

132 extinction coefficients were applied according to iFitzpatrick I (119991 ). with the largest source 

133 of error resulting from deredenning. A summary of the UVOT analysis results is given in 

134 Table 1. 



135 2.5. Optical MDM 

136 The region around RX J0648. 7+1516 was observed in the optical B, V, and R bands 

137 with the 1.3-m McGraw-Hill Telescope of the MDM Observatory on four nights during 2010 

138 April 1-5 (MJD 55287-55291). Exposure times ranged from 90 sec (R-band) to 120 sec 

139 (B-band). Each night, five sequences of exposures in B, V, and R were taken. The raw data 
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wo were bias subtracted and flat-field corrected using standard routines in Aperture 

141 photometry is performed using the IRAF package DAOPHOT on the object as well as five 

142 comparison stars in the same field of view. Calibrated magnitudes of the comparison stars 

143 are taken from the NOMAD catalog^, and the magnitudes of the object are determined using 

144 comparative photometry methods. For the construction of th e SED points, the magnitudes 

145 are extinction corrected based on the ISchlegel et al.l (119981 ) dust map with values taken 
we from NASA Extragalactic Database (NED)@ : A B = 0.618, A v = 0.475, and A R = 0.383. 
147 These data (summarized in Table 1) are used to constrain the modeling shown in this work, 
we although the same conclusions result with the UVOT points as model constraint. 



3. Spectroscopic Redshift Measurements 



150 
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Two spectra were obtained during the nights of UT 2010 March 18 and 2010 November 6 
(MJD 55245 and 55506, respectively) with the KAST double spectrograph on the Shane 3-m 
Telescope at UCO/Lick Observatory. During the first night, the instrument was configured 
with a 600/5000 grating and 1.5" long slit, covering 4300 — 7100 A. A single 1800 second 
exposure was acquired. During the night of November 6, another 1800 second exposure was 
acquired with a 600/4310 grism, D55 dichroic, a 600/7500 grating and 2" long slit, covering 
the interval 3500 — 8200 A. The data were reduced with the LowRedux pipelined and flux 
calibrated using a spectro-photometric star. The flux calibration is uncertain due to non- 
photometric conditions. Inspection of the March spectrum reveals Ca H+K absorption lines 
at redshift z = 0.179. This redshift is confirmed in the second spectrum at higher signal-to- 
noise (S/N) (S/N - 20 in the blue and S/N ~ 50 in the red) where Ca H+K, G band, Mgl 
AAA 5168,5174,5184 and Nal AAA 5891,5894,5897 absorption lines with equivalent width 
< 5 A are detected (see Figure [2] and Table 2 for details). No Ca H+K break is observed. 



Marcha et al. 


(1996 


) and 


Healev et al. 


(2007) 



4 http : //www.noao . edu/credit .html 
J http : //www.nof s .navy .mil/nomad .html 
e http : //nedwww. ipac . caltech . edu/ 



http : //www . ucolick . org/^xavier/LowRedux/ index . html 
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4. Broadband SED Modeling 
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The contemporaneous multiwavelength data are matched with archival radio data from 
NED and are shown in Figure 3. Since the radio data are not contemporaneous they are 
shown only for reference. The synchrotron peak appears at a frequency greater than 10 16 Hz, 
representing the first subclassification of RX J0648. 7+1516, specifically as an HBL. These 
data are used to test steady-state leptonic and lepto-hadronic jet models for the broadband 
blazar emission. The absorption of VHE gamma ra ys by the extragalactic background light 
(EBL) is accounted for thr ough application of the iGilmore et al. I (120091 ) EBL model; the 
model of iFinke et al.l ( 120101 ) provides comparable results. 



Leptonic models for blazar emission attribute the higher-energy peak in the SED to 
the inverse-Compton scattering of lower-energy photons off a population of non-thermal, 
relativistic electrons. These same electrons are responsible for the lower-energy synchrotron 
emission making up the first peak. The target photon field involved in the Compton upscat- 
tering can either be the synchrotron photons themselves, as in synchrotron self-Compton 
(SSC) models, or a photon field external to the jet in the case of external Compton (EC) 
models. 



We use the equilibrium SSC model of lBottcher fc Chiangl ( 120021 ) . as described in lAcciari et al 



( 120091 ) . In this model, the emission originates from a spherical blob of relativistic electrons 
with radius R. This blob is moving down the jet with a Lorentz factor T, corresponding to 
a jet speed of Pre. The jet is oriented such that the angle with respect to the line of sight is 
9 ohs , which results in a Doppler boosting with Doppler factor D = (T[l — (3 r cos O^])^ 1 . In 
order to minimize the number of free parameters, the modeling is completed with o b s = 1/r, 
for which Y — D. 

Within the model, electrons are injected with a power-law distribution at a rate Q{^) = 
Qo r )~ q between the low- and high-energy cut-offs, 7^2. The electron spectral index of 
q = 4.8 required for the models applied in this work might be the result of acceleration 
in an oblique shock. While standard shock acceleration in relativistic, parallel shocks is 
known to produce a canonical spectral index of ~2.2, oblique magnetic-fiel d configurations 
reduce the acceleration efficiency and lead to much steeper spectral indices (IMeli fc Quenby 
2003c ISironi &: Spitkovskyll201lf ). The radiation mechanisms considered lead to equilibrium 
between the particle injection, radiative cooling and particle escape. The particle escape is 
characterized with an efficiency factor rj, such that the escape timescale t esc = rjR/c, with 
77 = 100 for this work. This results in a particle distribution streaming along the jet with 
a power L e . Synchrotron emission results from the presence of a tangled magnetic field B, 
with a Poynting flux luminosity of Lb- The parameters L e and Lb allow the calculation of 
the equipartition parameter e^ e = L B /L e . 
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an The top panel in Figure 3 shows the SSC model for RX J0648. 7+1516, with parameters 

202 summarized in Table 3. The model is marginally in agreement with the data only through use 

203 of parameters well below equipartition. The Fermi-LAT contemporaneous 95% confidence 

204 level upper limits in the energy ranges 1-3 GeV and 3-10 GeV are just above and below the 

205 one-zone SSC model predictions. Additionally, these SSC model predictions are above the 

206 2.3 year Fermi-LAT spectrum by more than a factor of 2, although this spectrum is not 

207 contemporaneous with the other data. Variation of the model parameters within physically 
2ob reasonable values does not provide better agreement between model and data. Generally, 

209 HBLs are well characterized by one-zone SSC models and hence these observations might 

210 suggest the existence of one or more additional emission mechanisms that contribute to the 
2n higher-energy peak. 

212 An external-Compton model is also used to describe the data. The EC model is a 

213 leptonic one-zone jet model with two additional parameters beyond the SSC parameters, 

214 the thermal blackbody temperature Tec and radiation energy density uec of the external 

215 photon field, which is assumed to be isotropic and stationary in the blazar rest frame. The 

216 EC model provides a better representation of the SED, as can be seen in the middle panel 

217 of Figure 3, with the parameters listed in Table 3. 

218 A lepto-hadronic model is also applied to the data. Within this model, ultrarelativistic 

219 protons are the main source of the high-energy emission through proton synchrotron radi- 

220 ation and pion produ ction. The resulting spectra of the pion decay products are evaluated 



22i with the templates of iKelner fc Aharonianl (120081 ) . Additionally, a semi-analytical descrip- 



222 tion is used to account for electromagnetic cascades initiated by the internal 77 absorption 

223 of multi-TeV photons by both the 7 r° decay photon s and the synchrotron emission of ultra- 



relativistic leptons, as explained in iBottcherl (120101 ) . Similar to the particle populations in 



225 the leptonic models described above, this lepto-hadronic model assumes a power-law distri- 

226 bution of relativistic protons, 72(7) oc r y~ q between a low- and high-energy cut-off, £'™ m ' max . 

227 This population of relativistic protons is propagating along the blazar jet and has a total 

228 kinetic luminosity of L p . The lepto-hadronic modeling results are above €b p equipartition 

229 and are shown in the bottom panel of Figure 3 with parameters (including energy partition 

230 fractions e^ p = Lb/L p and e ep = L e /L p ) summarized in Table 3. 

231 In conclusion, multiwavelength followup of the VERITAS detection of 1FGL J0648. 7+1516 

232 has solidified its association with RX J0648. 7+1516, which is identified as a BL Lac object 

233 of the HBL subclass. Other contemporaneous SEDs of VHE-detected HBLs can be well 

234 described by one-zone SSC models close to equipartition, while for RX J0648. 7+1516 this 

235 model provides a poor representation with parameters below equipartition. The addition 

236 of an external photon field for Compton up-scattering in the leptonic paradigm provides 
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237 a better representation of the gamma-ray (Fermi and VERITAS) data. Alternatively, a 

238 lepto-hadronic model is successful in characterizing the higher-energy peak of the SED with 

239 synchrotron emission from protons. Both of these latter models require super-equipartition 

240 conditions. 
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1 10 10 2 Energy (GeV) 



Fig. 1. — Top: The differential photon spectrum of RX J0648. 7+1516 between 200 and 650 
GeV measured by VERITAS between 2010 4 March and 15 April (MJD 55259-55301). The 
solid line shows a power-law fit to the measured flux derived with four equally log-spaced 
bins and a final bin boundary at 650 GeV, above which there are few on-source photons. 
A 99% confidence upper limit evaluated between 650 GeV and 5 TeV assuming a photon 
index of 4.4 is also shown. The shaded region shows the systematic uncertainty of the fit, 
which is dominated by 20% uncertainty on the energy scale. Bottom: The differential photon 
spectrum of RXJ0648. 7+1516 as measured by Fermi-LAT over 2.3 years between 2008 Aug 
5 and 2010 Nov 17 (MJD 54683-55517, grey circles) with the highest energy bin containing a 
95% confidence upper limit. Fermi-LAT upper limits from the VERITAS observation period 
are also shown (MJD 55259-55301, grey triangles). 
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Fig. 2.— Spectrum of RX J0648.7+1516 showing the Ca H+K, G-band, Nal and Mgl 
spectral features indicating a redshift of z — 0.179. Since the G-band arises in stellar 
atmospheres, we interpret this as the redshift for the host galaxy and not an intervening 
absorber. The blazar was observed at Lick Observatory using the 3— m Shane Telescope on 
6 November 2010. 
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Fig. 3. — The SED models applied to the contemporaneous multiwavelength data of 
RXJ0648. 7+1516. Fermi-LAT data points are shown for 2.3 years of data along with up- 
per limits extracted from data limited to the VERITAS observation period. The models 
shown here are constrained by the MDM points; modeling constrained by the UVOT data 
produces similar results. The top panel shows the synchr otron emission (dotte d line), the 
self-Compton emission (dashed) and the EBL-corrected (IGilmore et al. 1 120091 ) total one- 
zone SSC model (solid). The middle panel shows the synchrotron emission (dotted line), 
the self-Compton emission (dashed line), the external-Compton (dash-dotted line) and the 
EBL-corrected total EC model (solid). The bottom panel shows the electron (and positron) 
synchrotron emission (dotted line), the proton synchrotron emission (dash-dotted) and the 
EBL-corrected total lepto-hadronic model (solid). 
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Table 1. Analysis summary of the optical MDM (B, V, R) and Swift-VVOT (U, UVM2, 

UVW2) data. 



Band 


Date 




vF v Error 




(MJD) 


(Jy Hz) 


(Jy Hz) 


B 


55287 


7.47X10 11 


3.4X10 10 


B 


55289 


7.64x10" 


3.8xl0 10 


g 




c; 7K v i nil 


9 7 v 1 n^o 


B 


55291 


7.59X10 11 


3.4xl0 10 


V 


55287 


5.77X10 11 


3.5xl0 10 


V 


55289 


5.74X10 11 


3.7xl0 10 


V 


55290 


2.92X10 11 


1.6xl0 10 


V 


55291 


6.00X10 11 


3.6xl0 10 


R 


55287 


5.99X10 11 


4.2xl0 10 


R 


55289 


5.51X10 11 


3.7xl0 10 


R 


55290 


2.03X10 11 


1.5xl0 10 


R 


55291 


5.99X10 11 


4.3xl0 10 


U 


55288 


4.542X10 11 


6.8xl0 9 


U 


55292 


4.253X10 11 


6.3 xlO 9 


u 


55300 


3.856X10 11 


6.1xl0 9 


u 


55304 


3.737X10 11 


5.5xl0 9 


UVM2 


55274 


5.987X10 11 


8.8x10 s 


UVW2 


55273 


5.066X10 11 


7.9x10 s 



Table 2. Analysis summary of the VER J0648+152 Lick Observatory Kast spectrum from 

2010 November 5 (MJD 55505) 



Ions 


Rest Wavelength 

(A) 


Ccntroid a 

(A) 


FWHM 

(A) 


Redshift b 
Absorbed 


Observed E. W. c 

(A) 


Notes 


Ca II (K) 


3934.79 


4639.07 


20.7 


0.1789 


2.60 ±0.21 




Ca II (H) 


3969.61 


4678.26 


16.4 


0.1785 


2.47±0.19 




G band 


4305.61 


5077.46 


17.5 


0.1792 


1.70±0.18 




Mgl 


5174.14 


6102.32 


22.1 


0.1793 


2.35±0.20 


[1] 


Na I 


5894.13 


6951.66 


23.0 


0.1794 


2.48±0.15 


[2] 



a Based on Gaussian fit 

b Measured from line centroid 

c Error is only statistical 



Note. — [1] Blandcd with Mg I 5168.74 Mg I 5185.04 [2] Blanded with Na I 5891.61 and Na I 
5897.57 
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Table 3. SED Modeling Parameters: Summary of the parameters describing the 
emission-zone properties for the SSC, EC and lepto-hadronic models. See text for 

parameter descriptions. 



Parameter 


SSC 


External Compton 


Lepto-Hadronic 


L e [erg s _1 ] 


7.5 x 10 43 


4.9 X 10 41 


4.9 X 10 41 


71 


6.7 x 10 4 


8.2 X 10 4 


9 x 10 3 


72 


10 6 


10 6 


5 x 10 4 


q 


4.8 


4.8 


4.8 


B[G] 


0.14 


0.1 


10 


r = d 


20 


20 


15 


Tec [K] 




10 3 




«BC [erg cm -3 ] 




7.0 x 10~ 8 




L p [erg s _1 ] 






4.9 x 10 41 


E™ in [GcV] 






10 3 


E™ ax [GeV] 






1.5 x 10 1() 


Qp 






2.0 


£Be 


0.16 


41 


1.7 x 10 4 


CBp 






4.2 


^ep 






2.5 x 10~ 4 


t™ n [hr] 


1.1 


10.9 


7.2 



